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ABSTRACT 

We are conducting COLD GASS, a legacy survey for molecular gas in nearby galaxies. 
Using the IRAM 30m telescope, we measure the CO (1-0) line in a sample of ~ 350 
nearby {Dl ^ 100 - 200 Mpc), massive galaxies (log(M*/M0) > 10.0). The sample is 
selected purely according to stellar mass, and therefore provides an unbiased view of 
molecular gas in these systems. By combining the IRAM data with SDSS photometry 
and spectroscopy, GALEX imaging and high-quality Arecibo HI data, we investigate 
the partition of condensed baryons between stars, atomic gas and molecular gas in 
0.1 — lOL* galaxies. In this paper, we present CO luminosities and molecular hydrogen 
masses for the first 222 galaxies. The overall CO detection rate is 54%, but our survey 
also uncovers the existence of sharp thresholds in galaxy structural parameters such 
as stellar mass surface density and concentration index, below which all galaxies have 
a measurable cold gas component but above which the detection rate of the CO line 
drops suddenly. The mean molecular gas fraction Mh2/M^ of the CO detections is 
0.066 ± 0.039, and this fraction does not depend on stellar mass, but is a strong 
function of NUV— r colour. Through stacking, we set a firm upper limit of Mh2/M^ = 
0.0016 ± 0.0005 for red galaxies with NUV— r> 5.0. The average molecular-to-atomic 
hydrogen ratio in present-day galaxies is 0.3, with significant scatter from one galaxy 
to the next. The existence of strong detection thresholds in both the HI and CO 
lines suggests that "quenching" processes have occurred in these systems. Intriguingly, 
atomic gas strongly dominates in the minority of galaxies with significant cold gas that 
lie above these thresholds. This suggests that some re-accretion of gas may still be 
possible following the quenching event. 

Key words: galaxies: fundamental parameters - galaxies: evolution - galaxies: ISM 
- radio lines: galaxies - surveys 
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1 INTRODUCTION 

Perhaps the most fascinating aspect of nearby galaxies is 
the intricately interwoven system of correlations between 
their global properties. These correlations form the basis of 
the so-called "scaling laws" , which are fundamental because 
they provide a quantitative means of characterizing how the 
physical properties of galaxies relate to each other. Galaxy 
scaling relations also provide the route to understanding the 
internal physics of galaxies, as well as their formation and 
evolutionary histories. 

We currently enjoy a diverse array of scaling laws 
that describe the stellar components of galaxies, for 
example the Tully- Fisher relation for spiral galaxies 
(|Tu11v ns"tolll97/tl. a nd th e fundamental plane for el- 
lipticals ('jorgen sen et al.|[l996l ). Both relations provide im- 
portant constraints on how these systems have assembled. 
However, few well-established scaling laws exist describing 
how the cold gas is correlated with the other global physical 
properties of galaxies. The only well- st udied relation is the 
Schmidt-Kennicutt star formation law (|Kennicuttlll998 '). re- 
lating the formation rate of new stars and the surface density 
of cold gas in disks. 

The reason why so few scaling laws involving cold gas 
and global galaxy properties, such as masses, sizes and 
bulge-to-disk ratios, exist in the literature, is the difficulty in 
acquiring suitable data. There are four general requirements 
on the data if the derived scaling laws are to be reliable: (1) 
homogeneous and accurate measurements of all the physical 
properties under consideration, (2) unbiased measurement 
of every property with respect to every other property, (3) 
sample selection that ensures large dynamic range of the 
various physical properties under consideration, (4) a large 
enough sample to define both the mean relation and the 
scatter about the mean. As we will describe, existing data 
sets do not, in general, meet all of these conditions. 

Line emission from the CO molecule was first de- 
tected in t he central parts and disks of nearby galaxies 35 
years ago (iRickard et al.1 1 19751 : ISolomon de Zafral Il975l : 
ICombes et al. 19771 ). A decade later, CO measurements ex- 
isted for ~ 100 galaxies (see the compilation of|Verter 19851), 
and in the following years, several larger systematic studies 
of molecular gas in nearby galaxies were performed. The 
largest effort was the FCRAO Extragalactic CO Survey 
("Young et allll995l ). which measured the CO J = 1 ^ 
line (hereafter, CO (1-0)) in 300 nearby galaxies. Because 
of a strong correla tion between infra r ed luminosity and 
CO l uminosity fe.g. ISanders et anil99ll : ISanders Mirabell 
Il996l \ a significant part of this early work was done tar- 
geting lumin ous infrared galaxies fe.g. 'Radford et 'al]|l99ll : 
ISolomon eF al. 1997). Even the FCRAO Survey, stih consid- 
ered as the reference for CO measurements in the nearby 
universe, targeted galaxies selected on infrared or B— band 
luminosity. 

Recognizing this bias toward "excepti onal" galax- 
ies (e .g. starbursts and interacting systems), iBraine et al.l 
(|l993l ) observed the CO 1^0 and 2^1 rotational tran- 
sition lines with the IRAM 30m telescope for a magnitude- 
limited sample of 81 normal spiral galaxies. Other attempts 
at meas uring molecular gas in n ormal galaxies include the 
work of iKennev Youna (| 19881 ) for Virgo cluster spirals. 



of Sagd (I1993II for nearby non-interacting spirals, and of 
iBoselh et al ] (|l997l ) for Coma cluster spirals. 

These pioneering studies constrained molecular gas 
properties in nearby galaxies as a function of mor- 
ihology (e.g. IWiklind Henkell Il989l : iThronson et ahl 
19891 ), star formation r a te or infrared lumino sity (e.g. 



Sanders k Mirabell Il985l: ICao Solomon! l2004l \ atomic 



conte nts (e.g. "Young & Knezek 1989), environment 
(e.g. Kenney k Young 1989; Casoh et al. 1991; Bosel li et al] 
I1997I ). and for resolved studies, position within galaxy 
disks. Highlights from these studies include the observations 
that molecular gas distributions decline monotonically with 
galaxy- centric radius unlike the atomic gas distributions, 
that IR- luminous galaxies are also CO-bright, with molec- 
ular gas concentrated within the inner kpc of these mostly 
interacting systems, and that the total gas mass fraction as 
well as the molecular-to-atomic ratio are functions of Hubble 
type. 

Nevertheless, most of the samples did not meet all of 
the criteria listed above that would allow for accurate scal- 
ing laws to be derived; some samples were biased towards a 
particular galaxy type (e.g. infrared-bright objects), some of 
the more unbiased samples did not cover enough parameter 
space (e.g. targeting only spiral galaxies), some samples suf- 
fered from aperture problems, some samples were too small, 
and attempts to combine different samples to remedy these 
problems led to inhomogeneous datasets (see also JT)). 

Recently, much effort has been put into obtaining ho- 
mogeneous and relatively deep high spatial resolution molec- 
ular gas maps covering the optical disks of nearby galax - 
ies (|Regan et al.ll200ll : lKuno et ahllioOTl : iLerov et al.ll2009l l. 
These samples are ex cellent for st udying star formation laws 
within galaxies (e.g. iBigiel et al . 2008), but the number of 
objects is too small to adequately define global scaling rela- 
tions. 

With reliable measurements of molecular gas for a large, 
unbiased sample of galaxies, it is possible not only to quan- 
tify scaling relations, but also to construct an accurate 
molecular gas mass function. Current es timates are based 
on inherently inhomogeneous samples (|Keres et al.l l2003l : 
lObreschkowH Rawlings 2009). We can also investigate the 
molecular gas properties of galaxy samples for which ded- 
icated surveys do exist, but where the number of objects 
studied has been very small, for example ear ly- type galax- 
ies (e.g. ICombes el al. "2007'; Krips et al."2010') , and galaxies 



with active nuclei (e. g. Heifer k Blitz 1993; Sakamoto et al.l 
I1999I : I Garcia- Burillo et al.l 120031 ) . A large unbiased sample 
of galaxies which can serve as a reference for such particular 
objects would also be very valuable. 

In this paper, we introduce COLD CASS, a new survey 
for molecular gas in nearby galaxies. Upon completion, it 
will have measured fluxes in the CO (1-0) line for a purely 
mass-selected sample of at least 350 galaxies. The sample 
contains galaxies with a wide range of Hubble types from 
star- forming spirals to "red and dead" ellipticals. With its 
new, large-bandwidth receivers, the IRAM 30m telescope is 
the instrument of choice to conduct a new large molecular 
gas survey, allowing the community to move from dedicated 
studies of particular types of galaxies, to larger systematic 
efforts. COLD CASS will provide a definitive, unbiased cen- 
sus of the partition of condensed baryons in the local Uni- 
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verse into stars, atomic and molecular gas in galaxies cover- 
ing over two orders of magnitude in luminosity. 

In 32]|4l we present an overview of the survey and of the 
sample selection, and describe the CO measurements and 
ancillary datasets. In JSl and [6] we present the first COLD 
GASS scaling relations, correlating molecular gas masses 
with global galaxy parameters including stellar mass and 
atomic gas mass. Throughout the paper, distance-dependent 
quantities are calculated for a standard flat ACDM cosmol- 
ogy with Hq — 70km s~^ Mpc~^, and we adopt a conversion 
factor from CO luminosity to H2 mass of aco = 3.2M0 (K 
km s~^ pc^)~^ (which does not account for the presence of 
Helium), unless otherwise specified. 



2 SURVEY DESCRIPTION AND SAMPLE 
SELECTION 

The conditions required to obtain reliable scaling laws and 
listed in JT] are routinely met by optically-selected samples 
of gala xies at low redshift. T he Sloan Digital Sky Survey 
fSDSS: Stmighton et al.ll2002l ) with its 5-band optical imag- 
ing campaign over a quarter of the sky and its follow-up 
spectroscopy of close to a million galaxies has facilitated the 
study of galaxy scaling relations studies at an unprecedented 
level of detail. 

At radio wavelengths, a series of large blind HI surveys 
have become possible thanks to a number of new multi-feed 
arrays. The most advanced of these, the Arec ibo Legacy 
Fast ALFA Survey fALFALFA: ICiov anelh et al.L 2QQ5). will 
have detected upon completion ^ 30, 000 galaxies out to dis- 
tances of 200Mpc. Although ALFALFA measurements are 
accurate, homogeneous, and unbiased, the survey is shallow, 
with the result that it does not probe a large dynamic range 
in Hl-to-stellar mass ratio for all but the very nearest galax- 
ies. For example, in the redshift range 0.025 < z < 0.05, 
the median value of Mhi /M^ for ALFALFA detections with 
M* > lO^^Mo is - 25%. 



2.1 GASS 

To overcom e this issue, the GA LEX Arecibo SDSS Sur- 
vey (GASS: ICatinella et al.ir201Ql l was designed to measure 
the neutral hydrogen content for a large, unbiased sample 
of ^ 1000 massive galaxies (M* > 10^° M©), via longer 
pointed observations. GASS is a large program currently 
under way at the Arecibo 305m telescope, and is producing 
some of the first unbiased atomic gas scaling relations in the 



some 01 the hrst unbiased atomic gas scalmg relations m the 
nearby universe (ICatinella et ahluOlOl : ISchiminovich et aP 
I2OIQI : iFabello et al.ll201Ql ~ 

Details about the GASS survey design, targe t selec- 
tion, a nd observing procedures are given in Catinella et al.l 
In short, the galaxies observed as part of GASS are 
selected at random out of a larger parent sample of galaxies 
that meet the following criteria: 

(i) They are located within the area of overlap of the 
SDSS spectroscopic survey, the ALFALFA survey, and the 
projected footprint of the GALEX Medium Imaging Survey 
(MIS). 

(ii) They lie in the redshift range 0.025 < z < 0.05. 
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Figure 1. Distribution of sources observed as of 25 Oct. 2010 
(filled blue histograms), compared to the proposed final COLD 
GASS sample (filled gray histograms). The solid line in each panel 
shows the distribution of objects in the superset of ~12000 galax- 
ies out of which the GASS sample is extracted, scaled down to 
the number of objects in the COLD GASS master list (350). 
The GASS sample is selected as to produce a flat distribution 
in log M* . 



(iii) They have a stellar mass in the range 10^*^ < 
M*/M0 < 10^^ ^ 

The GASS sample is selected out of this parent sample to 
produce a flat log M* distribution. No other selection criteria 
on colour, morphology, or spectral properties for example are 
applied. This sample therefore provides us with a complete 
picture of how the cold atomic gas relates to other properties 
such as stellar mass, luminosity, stellar surface mass density 

and colour. 

ICatineha et"all (l2010h present the first GASS data re- 
lease, which includes ^ 20% of the final sample. They show 
that there exist strong anti-correlations between the atomic 
gas mass fraction and stellar mass, stellar mass surface den- 
sity and NUV— r colour. GASS also aims at studying the 
galaxies that are transitioning between a blue, star-forming 
state and a red passive state (and vice- versa). These are 
identified as outliers from the mean scaling relations. The 
ultimate goal is to understand the physical processes that 
affect the gas content of these galaxies (e.g. accretion or 
quenching) and in turn the star formation process. 



2.2 COLD GASS 

We are in the process of constructing a CO Legacy Database 
for the GASS survey (COLD GASS), measuring the molecu- 
lar gas content of a significant subsample of the GASS galax- 
ies. We will then be able to quantify the link between atomic 
gas, molecular gas and stars in these systems. 

COLD GASS is designed to meet all the requirements 
to produce reliable gas scaling relations: 
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(i) Galaxies in our redshift range (0.025 < z < 0.05) have 
angular diameters that are small enough to enable accurate 
measurement of the total CO line flux with a single IRAM 
30m pointing for the majority (80%) of the galaxies. For 
the remaining objects, we recover the total flux by adding a 
single offset pointing (see ^4.41 for details). 

(ii) In the mass range we study (10^° < M*/M0 < 
10^^'^), the me tallicity of the galaxies is around solar 
( iTremonti et al ]f2004). The CO line flux therefore provides 
a reasonably accurate measurement of the total molecular 
gas content using a single conversion factor, aco- 

(iii) The ^ 350 targets are selected at random from the 
CASS survey, the sample is therefore unbiased. 

(iv) We integrate until the CO line is detected, or until we 
reach an upper limit in molecular gas mass to stellar mass 
ratio (fH2 = MhJM,) of - 1.5%. 

(v) Upon completion of the survey, the sample size of at 
least 350 galaxies will be large enough to determine accu- 
rately a set of scaling laws involving three parameters and 
to measure the scatter around these relations. 

Distributions of some basic parameters of the COLD 
CASS sample are shown in Figure [1] and compared to 
a purely volume-limited superset of galaxies. The COLD 
CASS targets are selected randomly from the CASS survey 
and therefore share the flat log M* distribution designed for 
that survey to ensure even sampling of the stellar mass pa- 
rameter space. As seen in Figured] a purely volume- limited 
sample is richer in low mass galaxies, and poorer in high 
mass systems. The uniform mass distribution also has the 
effect of flattening the colour distribution and reducing the 
number of galaxies with low stellar mass surface densities. 
We note that when deriving scaling relations and calculating 
sample ave rages, we statist i cally c orrect for this "mass bias" , 
as done bv lCatinella et all (|2010h for the CASS sample (see 

m ■ 



3 SDSS, GALEX AND ARECIBO 
OBSERVATIONS 

3.1 Optical and UV measurements 

Parameters such as redshifts, sizes, magnitudes, and Galac- 
tic extinction factor s are retrieve d from the database of 
SDSS DR7 ( Abazaii airet al.ll2009l l. The UV data are taken 
from the GALEX AW-sky a nd Medium Im aging surveys (AIS 
and MIS, respectively, see lMartin et al.t 2005). 

The SDSS and G ALEX images are reprocessed follow- 
ing 'Wang'et^l] in order to obtain accurate aperture 
photometry. The process includes registering the images and 
smoothing them to a common PSF. The SDSS r— band im- 
ages are convolved to the resolution of the UV imaging be- 
fore Sextractor is used to calculate magnitudes in consistent 
apertures, therefore ensuring that measurements in different 
bands represent similar physical regions of the galaxies. The 
derived NUV— r colours are corrected for Galact ic extinc- 
tion using the prescr iption of IWvder et al.l (|2007[ ) (see also 
ICatinella et al.ll2010h . 

Stellar masses are calculated from the SDSS photome - 
try using the SED-fitting technique of Salim et al.1 (|2QQ7|), 
assuming a Chabrier IMF (IChabr ier 2003). A variety of 
model SEDs from the iBruzual Charloti (|2003|) library are 



fitted to each galaxy, building a probability distribution for 
its stellar mass. The stellar mass assigned to a galaxy is then 
the mean of this distribution, while the measurement er- 
ror is estimated from its width. The systematic uncertainty 
between different technique to derive photometric stellar 
m asses from SDSS mea surements is < O.ldex, as estimated 
bv button eTaJ] (|2010l ). 

The main optical- and UV-derived parameters used 
throughout this paper are presented in Table [1] for all galax- 
ies within the present COLD CASS data release. Column 1 
and 2 give the CASS and SDSS ID numbers, respectively, 
column 3 gives the optical redshift from SDSS spectroscopy, 
while column 4 lists the stellar mass and column 5 the stellar 
mass surface density, which we calculate as: 



where i?5o,z is the 2;— band 50% flux intensity petrosian ra- 
dius, in kiloparsecs. In column 6, we give the ^— band opti- 
cal diameter {D25), and in column 7 the concentration in- 
dex (C = R90/R50, where R50 and R90 are from r— band 
photometry). Finally, columns 8 and 9 present the NUV— r 
colour and the r— band model magnitude, both corrected for 
Galactic extinction. 



3.2 HI masses 

Details of t h e H I observations are described in 
ICatinella et al.1 (|2010l ). so we only provide a brief overview 
here. The survey builds upon existi ng HI datab ases: the 
Cornell digital HI archive (Springo b et a lI I2005I ) and the 
ALFALFA survey (Giovan elh et al.1 l2005l ). HI data for 
about 20% of the GASS sample (the most gas-rich objects), 
can be found in either of these sources. For the rest of 
the sample, observations are carried out at the Arecibo 
Observatory. Integration times are set such as to detect 
HI gas mass fractions {/hi = Mhi/M^) of 1.5% or more. 
Observations are carried out using the L— band Wide 
receiver and the interim correlator, providing coverage of 
the full frequency interval of the GASS targets at a velocity 
resolution of 1.4 km s~^ before smoothing. Data reduction 
includes Hanning smoothing, bandpass subtraction, radio 
frequency interference (RFI) excision, flux calibration and 
weighted combination of individual spectra. Total Hl-line 
fluxes, velocity widths and recessional velocities are then 
mea sured using linea r fitt ing of the edges of th e HI profiles 
(e.g. ISpringob et al.ll2005l : ICatinella et ani2007l ). 



4 IRAM OBSERVATIONS 
4.1 Observing procedure 

Observations are carried out at the IRAM 30m telescope. 
We use the Eight Mixer Receiver (EMIR) to observe the 
CO(l-O) line (rest frequency, 115.271 GHz) . The CO 1 ^ 
transition traces well the entire molecular gas contents of 
th e galaxies at MH2) > 10^ cm~^ (see e.g. the Appendix 
in lTacconi et al.ll2008l ). In the 3mm band (E090), EMIR of- 
fers two sidebands with 8 GHz instantaneous bandwidth per 
sideband and per polarisation. With a single tuning of the 
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Table 1. Optical and UV parameters of the COLD GASS galaxies"^ 



GASS ID 


SDSS ID 


ZSDSS 


[log Mo] 


[logMokpc-^] 


1^25 

n 




NUV-r 
[mag] 


r 

[mag] 


11956 


J000820. 76+150921. 6 


0.0395 


10.09 


8.48 


22.5 


2.15 


3.04 


16.28 


12025 


J001934. 54+161215.0 


0.0366 


10.84 


9.13 


34.3 


3.03 


5.93 


14.73 


12002 


J002504. 00+145815. 2 


0.0367 


10.48 


9.41 


24.2 


3.17 


6.25 


15.46 


11989 


J002558. 89+135545. 8 


0.0419 


10.69 


9.18 


23.7 


3.02 


5.79 


15.13 


27167 


J003921. 66+142811. 5 


0.0380 


10.37 


9.14 


21.1 


2.77 


4.48 


15.49 


3189 


J004023. 48+143649. 4 


0.0384 


10.05 


7.92 


37.7 


1.96 


2.77 


15.65 


3261 


J005532. 61+154632. 9 


0.0375 


10.08 


8.57 


22.8 


2.54 


2.63 


15.48 


3318 


J010238. 29+151006. 6 


0.0397 


10.53 


8.98 


26.4 


3.05 


5.73 


15.21 


3439 


J010905.96+144520.8 


0.0386 


10.35 


8.78 


32.5 


2.90 


3.05 


15.48 


3465 


J011221. 82+150039.0 


0.0292 


10.19 


8.93 


28.7 


2.89 


3.63 


15.33 


3645 


J011501. 75+152448.6 


0.0307 


10.33 


8.93 


28.1 


2.71 


3.97 


15.11 


3509 


J011711. 65+132027.3 


0.0484 


10.81 


9.18 


31.1 


3.11 


4.14 


15.27 


3519 


J011728.11+144215.9 


0.0427 


10.74 


8.64 


34.2 


2.20 


3.68 


14.94 


3505 


J011746.76+131924.5 


0.0479 


10.21 


8.83 


17.7 


3.30 


4.92 


16.35 


3504 


J011823.44+133728.4 


0.0380 


10.16 


7.91 


37.7 


1.84 


2.85 


15.34 



^ Table 1 is published in its entirety in the electronic edition of the journal. A portion is shown here as example of 
its format and contents. 



receiver at a frequency of 111.4081 GHz, we are able to de- 
tect the redshifted CO (1-0) line for all the galaxies in our 
sample (0.025 < z < 0.05), within the 4 GHz bandwidth 
covered by the correlators. This single tuning procedure re- 
sults in enormous time savings of 15 minutes per source, 
and in an improved relative calibration accuracy. Also, the 
frequency range covered benefits from a considerably im- 
proved atmospheric transmission as compared to the C0(1- 

0) rest frequency. The second band is tuned to a frequency 
of 222.8118 GHz (E230 band), to cover the redshifted C0(2- 

1) line which falls within the available 4 GHz bandwidth for 
about 75% of our sample. We postpone the presentation and 
analysis of the CO (2-1) data until the survey is completed, 
to maximize the sample size. 

The wobbler-switching mode is used for all the obser- 
vations with a frequency of IHz and a throw of 180 " . The 
Wideband Line Multiple Autocorrelator (WILMA) is used 
as the backend, covering 4 GHz in each linear polarisation, 
for each band. WILMA gives a resolution of 2 MHz (~ 5 
km s~^ for the 3mm band). We also simultaneously record 
the data with the 4MHz Filterbank, as a backup. 

Observations for this first data release were conducted 
between 2009 December and 2010 October. Atmospheric 
conditions varied greatly, with an average of 6mm of pre- 
cipitable water vapor (PWV). We also fold into this catalog 
15 galaxies observed in 2009 June, as part of a pilot program 
designed to test the feasibility of the survey. These galax- 
ies were selected to be Hl-rich (Mhi/M^ > 0.1), but this 
selection bias does not affect the overall sample. 



4.2 Observing strategy 

Observations are carried out in fixed observing blocks and 
as poor- weather backups for higher frequency programs. We 
accommodate to the changing weather conditions by mak- 
ing "real time" decisions on targets. We observe the bluer 
galaxies (generally CO-luminous) under poorer weather con- 
ditions. These galaxies require on average an rms sensitivity 
of 1.7 mK per 20 km s~^-wide channel to achieve a reliable 
detection of the CO(l-O) line with S/N > 5. As seen in Fig- 
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Figure 2. Fraction of galaxies with a detection in the CO (1-0) 
line as a function of (a) stellar mass, (b) stellar mass surface den- 
sity, (c) concentration index (defined as Rqq/R^q, the ratio be- 
tween the r— band radii encompassing 90% and 50% of the light), 
and (d) NUV— r colour. The results are shown in equally popu- 
lated bins, each containing 37 galaxies. The gray shaded region 
shows the overall detection rate of 53.6%, down to 47.3% if all 
tentative detections are excluded. The downward error bars show 
the effect of excluding tentative detections in each individual bin. 
In panels (b-d), the vertical dotted line indicates the critical value 
where the detection rate suddenly drops below ~ 80%. 



ure[2jl, galaxies with colour NUV— r< 4.5 have a detection 
rate greater than 80%. When the atmospheric water vapor 
level is low, we preferentially observe the redder galaxies, 
which have a very low detection rate (Fig. [2li). In order to 
set firm upper limits for these galaxies, we require low noise 
values to reach our integration limit of Mh2/M^ = 1.5%. 
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Under good observing conditions, sensitivity to this mini- 
mum gas fraction, or an absolute minimum rms of l.lmK 
(per 20km s~^-wide channel), is reached within 1-1.5 hour. 
This absolute minimum rms is imposed in order to keep the 
integration time per galaxy < 2 hours, and translates in a 
detection limit that is higher than the nominal value of 1.5% 
for galaxies with log M^/M© < 10.6. 

The efficiency of the observations is also maximized by 
our single tuning approach (see ^4.ip . and by the fact that 
the galaxies are mostly concentrated in a declination strip 
with 0° < S < +15°, as shown in FigureO Not only does this 
allow us to move quickly from one source to the next with- 
out repeating pointing correction measurements, it makes it 
possible to almost always observe at elevations larger than 
45°, minimizing atmospheric opacity. 

4.3 Data reduction 

The data are reduced with the CLASS software. All scans 
are visually examined, and those with distorted baselines, in- 
creased noise due to poor atmospheric conditions, or anoma- 
lous features are discarded. The individual scans for a sin- 
gle galaxy are baseline-subtracted (first order fit) and then 
combined. This averaged spectrum is finally binned to a res- 
olution of ^ 20km s~^, and the standard deviation of the 
noise per such channel is recorded (arms)- 

Flux in the CO (1-0) line is measured by adding the 
signal within an appropriately defined windowing function. 
If the line is detected, the window is set by hand to match 
the observed line profile. If the CO line is undetected or very 
weak, the window is set either to the full width of the HI 
line {W50hi) or to a width of 300 km s~^ in case of an HI 
non-detection. For the non-detection, an upper limit for the 
flux of 5eobs (see Eq. [5]) is set. 

The central velocity and total width of the detected 
CO lines are then measured using a custom-made IDL in- 
teractive script. The peaks of the signal are identified, and 
a linear fit is applied to each side of the profile between 
the 20% and 80% peak flux level. The width of the line, 
WbOco , is then measured as the distance between the points 
on each of the fits corresponding to 50% of the peak inten- 
sity. The recession velocity is ta ken as the midp oint of this 
line. This metho d is described in lSpringob et al. l (2005) and 
ICatinella et al.l (|2007l V and is also used to measure the HI 
line- widths of the GASS sample. 



4.4 Aperture Corrections 

Because the galaxies targeted are at a distance of at least 
100 Mpc {z > 0.025), most of them can be observed with a 
single pointing of the IRAM 30m, which has a beam with 
a FWHM of 22^^ at a wavelength of 3mm. However some of 
the galaxies have optical diameters in excess of this, and an 
aperture correction needs to be applied. 

We derive aperture correct ions usin^ a set of nearby 
galaxies with accurate CO maps (|Kuno et al.ll2007l ). We sim- 
ulate the impact of observing galaxies with the IRAM beam 
by taking each of the maps, placing it at different redshifts 
in the range 0.025 < z < 0.05, and computing the ratio 
between the flux as would be measured by a 22^^ Gaussian 
beam to the total flux in the map {Sco.obs/ Stot)- We find 




20 40 60 80 100 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
D25 [arcsec] S„ff/Sco.obs 



Figure 4. Aperture corrections for the IRAM 30m observations 
of the COLD GASS sample, (a): Ratio between the flux recovered 
by a central pointing with the IRAM 22^^ beam to the total flux 
of the galaxy (Sco.obs/ ^tot), as a function of optical diameter 
(^25)- Small gray points indicate the results from the simulated 
observations of a samp le of 40 near by spiral galaxies with high 
quality CO (1-0) maps (|Kuno et al. l 2007) when they are placed 
at various redshifts in the range 0.01-0.05 (each set of points con- 
nected by a broken line corresponds to one of the galaxies, placed 
at different redshifts). A single central pointing recovers at least 
60% of the flux in galaxies with D25 < 40^^ The histogram indi- 
cates the size distribution of the COLD GASS sample (with the 
distribution of blue galaxies, NUV— r<4.3, indicated in blue), (b): 
With an additional offset observation 0.75 beam (i.e. 16^^) from 
the centre of the galaxy along the major axis (Soff), it is possible 
to derive the fraction of CO (1-0) line flux detected at the center 
and from that estimate the total molecular gas mass in galax- 
ies with D25 > 40^^ The large red symbols in panel (a) indicate 
the aperture corrections estimated using this method for the 25 
COLD GASS galaxies for which we performed offset pointings to 
date, and the dashed line is the size threshold (40^^) for a galaxy 
to require an offset pointing. 

that a single central observation recovers most (> 60%) of 
the CO line flux in galaxies with D25 < 40^^ Results are 
shown in Figure |4^. Based on the best fit to these data, we 
apply the following aperture correction: 

Sco,cor = Sco,obs/{l.OM - OmSD25 + 2.0 X 10-'i:>^5) (2) 

where Sco.obs is the observed flux in the central pointing 
and S CO, cor the extrapolated total flux. 

For galaxies with optical diameters larger than 40^^ 
however, there is a significant scatter in the Sco,obs/ Stot 
ratio, and additional information is required to recover the 
total CO line flux. In the right panel of Figure 2] we show 
that with a single offset pointing at three-quarter beam from 
the central position (0.75x22^^) along the major axis, the 
total flux can be recovered with much better accuracjQ. An 
offset pointing at a full beam also does well, but the mean ra- 
tio Sof f / Sco,obs then drops from 33% to 15%. We adopted 
the three-quarter beam offset as a compromise between the 
requirements for independent flux measures and a modest 
fraction of our total observing time going into off-center 
pointings. 



^ During the pilot observations of 2009 June, we took offset point- 
ings one full beam from the central position. These galaxies are 
identified clearly by asterisks in Col. 6 of Tab. O and aperture 
corrections are performed using a version of Eq. [3] appropriate for 
these larger offsets. 
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Figure 3. Sky distribution of the COLD GASS sample. Galaxies are selected in the area of intersection of the SDSS, GALEX 
MIS and ALFALFA HI surveys. Gray dots show all 12000 galaxies satisfying the GASS selection criteria (see within the 

SDSS+GALEX footprints, and the circles represent the 222 COLD GASS galaxies that are part of this data release. 



Our requirement to perform an offset pointing is that 
a galaxy has (1) a large angular size {D25 > 40") and (2) a 
bright CO line in the central pointing, such that a detection 
in the offset pointing can be made with an hour of integra- 
tion or less. For these galaxies, the flux in the central point- 
ing is corrected based on the ratio between the flux in the 
offset pointing to that in the center {foff = Sof f / Sco.obs), 
based on the best fit to data in Figure HJ): 



Sco,^ 



^co,ob./(1.587 - 3.36I/0// + 2.107/0' 



ff) 



(3) 



If a galaxy is larger than 4Q" but the CO line in the central 
pointing is weak or undetected, then the measured central 
flux is corrected according to Eq. [2] 

So far, we have performed offset pointings for 25galax- 
ies, that met the requirements listed above (see Figs. IA14I 
and I A 15")) . We used Eq. [3] and Fig. UJ) to infer the fraction 
of the total flux that was measured in the central pointing 
{Sco.obs/ Stot), and we show in Fig. how this ratio de- 
pends on D25 and how it compares to the range of values 
predicted by our simulated observations. 

4.5 Mh2 and associated error budget 

After correcting the CO (1-0) line fluxes for aperture ef- 
fects using either Eq. [2] or Eq. [S] whichever case applies to 
each galaxy, we compute the total CO luminosities following 
ISolomon et aD (|l997l ): 



L'co = 3.25 X 10^ Sco,coriy;j,lDl{l + z)-\ 



(4) 



where Sco cor in units of Jy km s ^is the integrated line 
Vobs is the observed frequency of the CO (1-0) line in 
GHz, Dl is the luminosity distance in units of Mpc, and 
L'co is the CO luminosity in [K km s~^ pc^]. 

The total molecular hydrogen masses are then calcu- 
lated as Mh2 = L^Qaco- We adopt a constant Galactic 
conversion factor of aco = 3.2 M0(K km s~^ pc^)~^, which 
does not include a correction for the presence of Helium. Our 
choice of aco is roughly the mean of values estimated in the 
Milky Way and in nearby galaxies (e.g. IStrong &; Mattoxl 

^ Calculated from antenna temperature units using the conver- 
sion S/T* = 6.0Jy/K, specific for the IRAM 30m at our observing 
frequency of lllGHz. 



19961: iDame et al 11200 lUBlitz et al.l[2007l : iDraine et all2007l : 
Hever et al.ll2009l : lAbdo et al.ll201o[ l. The virial method used 



to measure aco has been validated by other independent 
techniques such as 7— ray observations, and shown to also 
hold for the ensemble average of the virialized clouds of en- 
tire galaxies, as long as the factor n(H2)°'^/T is constant 
througho ut the galaxy and the CO line is optically thick 
(|Dickman et al.lll986l : iTacconi et al.ll2010h . 

This constant value for aco has been shown to hold for 
galaxies in the Local Group, when the metallicit y goes from 
solar down to SMC values (Bola tto et alll2008l ). It may be 
that the conversion factor is instead a functio n of a pa rame- 
ter such as gas surface density or met allicity (|Tacconi et al.l 
I2OO8I : lObreschkow k Rawhngsl l2009l ). Based on the well- 
known metallicity-luminosity relation, iBoselh et al.] (|2002l ) 
for example proposed a luminosity-dependent conversion 
factor. However, such a prescription has yet to be obser- 
vationally or theoretically validated. Furthermore, it has 
since been shown t hat the mass-metallic ity relation is even 
more fundamental (iTremonti et al.ll2004[ ). At stellar masses 
above ^ lO^^'^M©, the relation flattens out, and therefore 
the metallicity of our COLD GASS galaxies is expected to 
be ^solar with little variations across the sample. This ex- 
pectation i s confirmed by our long- slit spectroscop y mea- 
surements (|Moran et al.ll2010l : ISaintonge et al.ll201lh . Based 
on this fact, adopting a constant Galactic conversion fac- 
tor aco — 3.2 M0(K km s~^ pc^)~^ is the simplest yet 
most justified assumption we can make at this point. We 
however consider the variat ions in aco by a factor of ^2 
at fixed metallicity (Bolatt o et ahllioosh as the systematic 
uncertainty on our values of • 

We calculate the formal measurement error on the ob- 
served line flux, Sco.obs, as: 



^obs 



sW50c 



y^W50coA% 



(5) 



where arms is the rms noise per spectral channel of width 
Awch = 21.57km s-\ and W50co is the FWHM of the 
CO (1-0) line. This definition takes into account that for a 
given total flux, the S/N per channel is largest when the 
emission line is narrower. Considering all detections, the 
mean fractional error {cobs/ Sco,obs) is 11%. 
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Other contributions to the error budget on Mh2 include 
a flux calibration error (10% at a wavelength of 3mm, un- 
der average atmospheric conditions) and the uncertainty on 
the aperture correction (which we estimate to be 15% based 
on Fig|4]for the median galaxy in our survey). The average 
pointing rms error is Apoint = 2^^, and also contributes to 
the uncertainty on the measured flux. Using the resolved 
CO (1-0) maps of iKuno et aP (|2007l ). we simulate the obser- 
vation of COLD GASS galaxies both with a beam perfectly 
centered on the object, and with a beam offset by Apoint 
from the center. We find that a 2" positional error gener- 
ates a 2.1% uncertainty on the measured line flux. Given 
that the redshift error is negligible compared to these other 
contributions, the fractional error on Lqq is obtained by 
adding in quadrature all error contributions to Sco,cor- We 
therefore find a mean error of < clco >— 20%. The random 
measurement error on our quoted values of Mh2 is then 
~ 20%, to which we add the systematic error coming from 
the uncertainty in the value of aco bringing the total error 
on logMi/2 to 0.3 dex. 



Column (2): rms noise per channel in mK, after bin- 
ning the spectra to a resolution of Awch — 21 km s~^. 

Column (3): Signal-to- noise ratio of the detected CO 
line, calculated as Sco.obs/^obs, where Cobs is calculated 
according to Eq. [5] In our analysis, we consider sources 
with 5'/A^ > 5 as secure, and the detections with S/N < 5 
as tentative. 

Column (4): Integrated CO line flux in Jy km s~^. The 
measured flux in antenna temperature units is converted 
to these units using the ratio S/T* = 6.0 Jy/K for the 
IRAM 30m telescope at our observing frequency of lllGHz. 

Column (5): Corrected total line flux in Jy km s~^, 
computed following Equations [2] and (3] 

Column (6): When available, the measured ratio be- 
tween the flux in the offset pointing to that in the central 
pointing, as described in ^4.41 



4.6 Stacking 

To extract information from the CO non-detections, we also 
perform a stacking analysis. We stack the spectra after con- 
verting them in "gas fraction" units following IPabello et al.l 
(|20ld l. who provide an extensive discussion of the merits of 
this technique. We convert each spectrum S into a scaled 
spectrum Sgf. 



SDi 



M*(l + z) 



(6) 



where the factor 0^^(1-1- z)~^ is a consequence of Eq.[31 This 
gives us scaled spectra in units of [mJy Mpc^ M©"^]. 

For any subsample of N galaxies with CO non- 
detection, we calculate a mean stacked spectrum: 



Ss 



EN ' 



(7) 



where w is the weight given to each galaxy in the stack, 
based on its stellar mass, in order to compensate for the fact 
that the log M* distribution of the COLD GASS sample is 
by design flatter than in a complete volume-limited sample 
(Figure [I^). 

To obtain < >, the mean molecular gas mass 

fraction for the family of N galaxies stacked, we measure 
a line flux or set an upper limit as described in ^4.3[ and 
multiply by the constant factors 3.25 x 10^ iy~j^^aco to obtain 
the dimensionless mass fraction. 



4.7 Catalog presentation 

In table H we summarize the results of the IRAM obser- 
vations of the first 222 COLD GASS galaxies (119 CO 
detections and 103 non-detections). The contents of the 
catalog are as follows: 



Column (7): Total molecular hydrogen gas mass. We 
adopt a constant Galactic conversion factor of aco = 3.2 
M0(K km s~^ pc^)~^, which does not include a correction 
for the presence of Helium. The numbers quoted are either 
measured masses in the case of detections, or 5cr upper 
limits for the non-detections (see Column 9). 

Column (8): Molecular gas mass fraction, fH2 = Mh2/M^^ 
and upper limits for non-detections. 

Column (9): CO emission line flag, set to 1 for detec- 
tions and 2 for non-detections. 



5 MOLECULAR GAS FRACTION SCALING 
RELATIONS 

5.1 Molecular gas fraction and global galaxy 
parameters 

We first investigate how the molecular gas mass fractiorjfl 
Uh2 = M//2/M*) varies across the COLD GASS sample. 
The values of fH2 range from 0.009 (G7286), at the very 
edge of the survey's detection limit, up to 0.20 for the most 
i72— rich galaxy in the sample observed so far (G41969). 
Overall, the molecular gas contents are more modest, aver- 
aging at < fH2 >= 0.066 =b 0.039 for the CO detections, 
and at 0.043 ± 0.022 when non-detections are included in 
the mean at the value of their upper limits. These sample 
averages are measured by weighting each galaxy differently 
according to its stellar mass, in order to compensate for the 
flat logM* distribution (see Section [5.21 for more details). 

In Figure O we investigate how the molecular gas mass 
fraction depends on stellar mass, stellar mass surface den- 
sity, concentration index and NUV— r colour. For all four 
relations, we quantify the strength of the dependence of fH2 
on the X-axis parameter in two ways. Firstly, we compute for 



Column (1): GASS ID. Galaxies are in the same or- 
der as in Table [1] to ease cross-referencing. 



^ We warn the reader that our definition of fH2 differs from 
some previous studies, where fH2 is also defined as Mh2/ ^dyn 
or Mh2/(M^ + Mhi + Mh2)- 
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Table 2. Molecular gas masses and CO (1-0) parameters for the COLD GASS galaxies" 



GASS ID 


(7 


S/N 


ScO,ohs 


ScO,cor 


foff 


Mh2 


Mh2/M^ 


Flag 




[mK] 




[Jy km s"-*-] 


[Jy km s"-*-] 




[log Mo] 






11956 


1.07 


2.35 


1.16 


1.25 




8.46 


0.024 


1 


12025 


1.06 










8.78 


0.009 


2 


12002 


1.18 










8.79 


0.021 


2 


11989 


1.07 










8.79 


0.013 


2 


27167 


1.17 










8.74 


0.023 


2 


3189 


1.24 


6.69 


3.19 


3.88 




8.93 


0.076 


1 


3261 


1.96 


8.57 


4.27 


4.62 




8.98 


0.080 


1 


3318 


1.03 










8.81 


0.019 


2 


3439 


1.03 










8.79 


0.028 


2 


3465 


1.17 


4.28 


2.89 


3.28 




8.62 


0.027 


1 


3645 


1.08 










8.62 


0.020 


2 


3509 


1.15 


7.24 


5.02 


5.80 




9.30 


0.031 


1 


3519 


1.53 


7.72 


5.51 


6.52 




9.24 


0.032 


1 


3505 


1.30 










8.93 


0.052 


2 


3504 


1.53 


9.19 


2.83 


3.44 




8.87 


0.051 


1 



" Table 1 is published in its entirety in the electronic edition of the journal. A portion is shown here 
as example of its format and contents. 



all secure detections the Pearson correlation coefficient of the 
relation, r. Secondly we perform an ordinary least squares 
linear regression of fH2 on the x-axis parameter for the se- 
cure detections, and report the standard deviation of the 
residuals about this best fitting relation, a, in logMi/2/M* 
units. 

There is a weak dependency of fH2 on stellar mass, M*, 
as seen in Fig. [5^, with a correlation coefficient r = —0.40, 
and a scatter of a = 0.292 dex, but some of this correlation is 
caused by the detection limit. The strength of the correlation 
of fH2 on stellar mass surface density, /i*, is similar with 
r = —0.33 and a = 0.301 dex (Fig. \Sjp>). There is however 
a striking difference in the behaviour of the non-detections; 
while the detection rate is roughly constant as a function 
of M*(^ 50%, see Fig. [2^), it is a strong function of /x*. 
There is a surface density threshold of /x* ^ 10^"^ M© kpc~^, 
below which the detection rate is almost 100% and above 
which it quickl y drop s. A similar threshold is observed by 
'Catinella et al.l (|2010l ) in the relation between the atomic 
gas mass fraction {fni = Mhi/M^) and /i* (see 31] for more 
details) . 

To test the impact of galaxy morphology on the molec- 
ular mass fraction, we plot in Fig. [St fH2 as a function 
of the concentration index, defined as C = R90/R50 where 
R90 and R50 are the radii enclosing 90% and 50% of the 
r— band flux, respectively. The concentration index is a 
good proxy from the bulge-to-total ratio, as would be re- 
covered by fu ll two-dimentional bulge/disk decompositions 
(jWeinmann et alj l2009l V The galaxies with a CO line de- 
tection show no correlation, with a correlation coefficient of 
—0.14, and a scatter about the best-fit relation (0.317 dex) 
as large as the scatter in fH2 itself (also 0.317 dex), suggest- 
ing that the presence of a bulge does not affect strongly how 
much of the mass in these galaxies is in the form o f molecular 
gas. A similar result was found bv lCatinella et all (|201fl) and 
iFabello et all (|201Q| 1 for the HI contents of similar galaxies. 
The only effect seen is again in the detection fraction, with 
the success rate of measuring the CO line dropping sharply 
from 100% for the most disk-dominated systems, to 0% for 
the most bulge-dominated ones. 
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Figure 6. Stacked spectrum of all galaxies with NUV— r> 5.0. 
Individually, the galaxies included are all non-detections in CO, 
and even their stacked average does not yield a detection, set- 
ting instead a stringent upper limit of < Mh2/M* >= 0.0016 ± 
0.0005. The vertical dotted line show an expected line width of 
~ 500km s"-*^ for such a stacked signal. 



The only parameter upon which the measured molecu- 
lar mass fractions does depend significantly is colour (r = 
-0.57, a = 0.260 dex, see Fig.[5]i). Since NUV-r colour is 
a proxy for specific star formation rate, this dependency of 
fH2 is expected, because star formation and molecular gas 
are known to be strongly correlated. It is interesting to note 
that not a single red sequence galaxy (NUV-r > 5.0) has 
a measurable molecular gas component {fH2 < 0.015 in all 
cases). To test this further, we stack all non-detections with 
NUV— r> 5 using the technique described in ^4.61 The re- 
sult is shown in Figure [G] Even in the stack, the red galaxies 
lead to a non-detection of the CO line, thus setting an even 
more restrictive upper limit on the molecular gas mass frac- 
tion of 0.0016 =b 0.0005, which at their median stellar mass 
corresponds to an upper limit on the molecular gas mass of 
^ 1.1 X 10^ M0. Molecular gas has been det ected in early- 
type galaxies (e.g. in the SAURON sample, ICombes et alj 
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Figure 5. Molecular gas mass fraction {fH2 = Mh2/M^) as a function of (a) stellar mass, (h) stellar mass surface density, (c) 
concentration index Rqq/R^o, a,nd(d) NUV— r colour . Detections in CO are shown as circles (filled: S/N > 5, open: S/N < 5), 
and upper limits obtained for non-detections as arrows. In panel (a), the dotted line represents the detection limit of the survey, 
which is a function of M* at the low mass end due to the hard limit of l.lmK for the observed rms (see Section [4. 21 for details). 
In each panel, we give the value of the Pearson correlation coefficient of the relation (r), and the standard deviation a about 
the best fit linear relation to all CO detections (dashed line). 



but at levels even below this limit. Our results are 
therefore not in contradiction with these studies. 



5.2 Mean scaling relations 

After having looked separately at the properties of CO de- 
tections and non-detections as a function of several global 
galaxy parameters (Fig. O, we combine all measurements 
into mean molecular gas mass fraction scaling relations. 

As described in ^2.21 the COLD GASS sample has been 
generated to have a stellar mass distribution that is flatter 
than it is in a purely volume- limited sample (see Fig. [1]). 
When building the scaling relations, we correct for this by 
wei ghing each poin t according to its stellar mass. Follow- 
ing jCannell^^et^l] (|2010l l, the galaxies are placed in bins 



of stellar mass of width 0.2 dex and assigned as a weight 
the ratio between the total number of galaxies in the unbi- 
ased volume-limited parent sample and the total number of 
COLD GASS galaxies within that same mass bin. In other 
words, low mass galaxies are given a higher weight in the 
computation of the mean scaling relations, because these 
galaxies are underrepresented in the COLD GASS sample 
compared to a volume-limited sample. 

Each scaling relation is computed and plotted using 
three different subsamples: (1) only the galaxies with CO 
detections are considered, (2) both the detections and non- 
detections are used, and the upper limit on fH2 is used for 
the non-detections, and (3) both the detections and non- 
detections are used, but this time a value of fH2 = 0.0 is 
assigned to the non-detections. In all cases, the weighted 
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Figure 7. Molecular gas mass fraction scaling relations. The weighted mean values of log fH2 in equally populated bins are 
plotted as a function of (a) stellar mass, (h) stellar mass surface density, (c) concentration index Rqo/R^q, a,nd(d) NUV— r 
colour. The blue circles include only CO detections. The red circles and green squares include both detections and non-detections, 
but differ in the way the latter are treated. We assign to the non-detections either the value set by the upper limit (red circles), 
or simply set them to zero (green squares). The correlation coefficients (r) and scatter around the best-fit relation are given for 
detections only (lower left corners) and for the "non-detections as upper limits" case (lower right corners). The gray symbols 
show the data for all COLD GASS galaxies within this data release (as presented in Fig. [5]). 



Table 3. Mean molecular gas fraction scaling relations' 



(1) Detections only (2) All, non-detections=limit (3) All, non-detections=0 



X 


xo 


m 


b 


m 


b 


m 


b 


logM* 


10.70 


-0.202 ±0.040 


-1.300 ±0.592 


-0.455 ± 0.069 


-1.607 ± 1.034 


-0.346 ± 0.203 


-1.552 ±3.038 


log/i* 


8.70 


-0.283 ±0.019 


-1.265 ±0.237 


-0.595 ±0.055 


-1.441 ±0.682 


-0.816 ±0.127 


-1.454 ± 1.587 


Rqo/R^o 


2.50 


-0.078 ± 0.048 


-1.258 ±0.176 


-0.447 ±0.049 


-1.420 ±0.185 


-0.644 ±0.145 


-1.388 ±0.529 


NUV-r 


3.50 


-0.197 ±0.012 


-1.275 ±0.066 


-0.219 ±0.013 


-1.363 ±0.077 


-0.293 ±0.013 


-1.349 ±0.069 



" The relations are parametrized as log t(^ep(H2)[yr ^] = Tn{x — xq) ± b. 
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mean of log fH2 is then calculated in equally populated bins 
of either M*, /i*, concentration index or NUV— r colour, 
and the relations are fitted, weighting galaxies according to 
their stellar mass. The resulting mean relations are plotted 
in Figure [71 with the error bars representing the uncertainty 
on the position of < log fH2 > in equally populated bins, 
as determined by bootstrapping: the error is the standard 
deviation in the value of < log fH2 > for 1000 resamples of 
the original data in each bin, with repetitions. The best-fit 
linear relations are also plotted and summarized in Table [S] 
Note that in case (3), because of the null values of fH2 for 
the non-detections, we cannot directly average and fit the 
values of \ogfH2 as we do for cases (1) and (2). Instead, we 
measure the logarithm of the mean value of fH2 in each bin, 
and fit these average values. This is the relation plotted in 
Figure [3 and given in Table [3] for case (3). 

The choice between setting the non-detections to their 
upper limits or to a constant value of zero only significantly 
affects the scaling relations at large values of /x*, C and 
NUV— r where few galaxies have a detected CO line. Both 
sets of scaling relations are included in Table [3l labeled as 
< log Mh2/Mq >iim and < \ogMH2/MQ >o for the upper 
limit and zero value cases, respectively. 

The mean molecular gas mass fraction is a roughly con- 
stant function of stellar mass, both for detections alone and 
when including non-detections, as shown in Fig. [7^. This is 
a consequence of the flat detection rate of the CO line as a 
function of M* . Adding non-detections turns the mostly- flat 
relations between fH2 and both /i* and C into monotoni- 
cally decreasing functions. The strongest correlation is still 
with NUV— r colour, both before and after including non- 
detections. Because NUV— r colour is a proxy for specific 
star formation rate, its correlation with fH2 is not surprising, 
given e.g. the Kennicutt-Schmidt relation. The correlation of 
fH2 with NUV— r could therefore be seen as a consequence 
rather than a cause. 

For the other parameters describing the underlying 
properties of the galaxies, independently of the current star 
formation rate (M*, /i*, concentration index), there is very 
little dependence of the measured values of ///2- In ^5.11 we 
showed that no matter their stellar mass, about half of our 
sample has detectable CO line emission, which converts into 
a molecular mass fraction of ^ 6% across the stellar mass 
interval sampled. On the other hand, while the molecular 
mass fraction of detected galaxies is mostly independent of 
/i* , the detection fraction is a strong function of that quan- 
tity. These observations trace a picture where: (1) the condi- 
tions required for the formation/consumption of molecular 
gas are a strong function of /x* (and concentration index) 
but not of M*, but (2) when these conditions are met, a 
roughly constant fraction of the stellar mass is found in the 
form of molecular gas. 



6 THE RELATIONSHIP BETWEEN HI AND 

Under the assumption that molecular gas forms out of lower 
density clouds of atomic gas, one might naively expect a 
tight correlation between Mhi and Mh2- The actual situa- 
tion is, however, quite different, as seen in Figure El Within 
the subsample of galaxies detected both in HI and CO, 



10,5 



10.0 



^ 9.5 
^ 9.0 
8.5 



8.0 



8.5 9.0 9.5 10.0 10.5 
log Mhi [Mq] 



Figure 8. Comparison between atomic and molecular hydro- 
gen gas masses. COLD GASS galaxies with detections of both 
the HI and CO lines are plotted as filled and open blue circles, 
for secure and tentative CO detections, respectively. The arrows 
shows limits in the cases of non-detecion of either or both the 
HI and CO lines. The best bisector linear fit to the detections is 
show as a solid line, and the ±1.5cr region around this fit marked 
with dashed lines. For comparison, we overplot as filled squares 
the integrated measureme nts for the HERAC LES galaxies with 
M* > IO^OMo, taken from lLerov et al.l (|2008l '). 



the fraction Mh2/Mhi varies greatly, from 0.037 (G13775) 
up to 4.09 (G38462). For galaxies with both CO and HI- 
detections, the correlation coefficient between log Mhi and 
log Mh2 is r = 0.37, indicating that the two quantities are 
only weakly correlated. A bisector linear fit to the same sub- 
sample reveals that on average, Mh2 is 0.295 times the value 
of Mhi, with a large scatter of 0.41 dex in \og{MH2/MHi) 
(see Fig.E)). 

To further investigate the relationship between Mh2 
and Mhi, we look into how the ratio between these two 
quantities varies as a function of different physical parame- 
ters. We define the molecular fraction as the ratio between 
the molecular hydrogen gas mass and the atomic gas mass 
of the system: 

In Figure [9l we plot Rmoi as a function of the same four 
global parameters as in the previous figures. Galaxies with 
non-detections in both HI and H2 are not plotted, as Rmoi 
is completely unconstrained in these cases, he best fit lin- 
ear relations, taking into account the weights correcting for 
the fiat logM* distribution, are measured in two ways: (1) 
including only the secure detections in HI and H2, and 
(2) including the non-detections in either of these quanti- 
ties as lower or upper limits, respectively. In Figure [9l we 
also report the Pearson correlation coefficients and the scat- 
ter around these relations. The best-fit scaling relations are 
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Figure 9. Molecular to atomic gas mass ratio as a function of (a) stellar mass, (h) stellar mass surface density, (c) concentration 
index Rqo/R^q, and(d) NUV— r colour . Detections in CO are shown as circles (filled: S/N > 5, open: S/N < 5), non-detections 
as upper limits. Non-detections in HI but with a CO detection are shown as lower limits. The mean scaling relations are also 
overplotted, with symbols as described in Fig. [Jl 



Table 4. Molecular-to- atomic gas mass ratio relations' 



(1) Detections only (2) All, non-detections=limit 



X 


xo 


m 


b 


m 


b 


logM* 


10.70 


0.425 ± 0.097 


-0.387 ± 1.464 


0.303 ± 0.084 


-0.489 ± 1.253 


log/i* 


8.70 


0.533 ±0.128 


-0.451 ± 1.562 


0.346 ±0.142 


-0.562 ± 1.772 


Rqo/R^o 


2.50 


0.152 ±0.031 


-0.475 ±0.113 


-0.106 ±0.098 


-0.571 ±0.367 


NUV-r 


3.50 


0.206 ±0.057 


-0.436 ±0.289 


0.075 ±0.032 


-0.575 ±0.169 



The relations are parametrized as log Mjj^ / ^Hl = m{x — xq) ± b. 



given the form Rmoi = m(x — xo) + with the parameters 
presented in Tabled 

The first panel shows that Rmoi is a very weakly increas- 
ing function of M* (r = 0.23). Both with and without the 
non-detections included, the slope of the relation is positive 



at the 3cr significance level (see Table [4|). As demonstrated 
previously, fH2 does not appear to depend significantly on 
M^(Fig.[5 tt), while .fa/ is a f airly strongly declining function 
of (^Catinella et al.ll201Ql V The combination of these two 
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trends produces the weak increase in Rmoi as a function of 
stellar mass. 



ICatinella et^ (|2010l l also reported a strong anti- 
correlation of fni with the stellar mass surface density, /x*. 
For galaxies with CO detections, the values of fH2 show 
a similar, but considerably weaker trend as a function of 
/i* . In their resolved study, Lero v et alj ([2008") show that 
the molecular-to-atomic ratio is a strong function of local 
properties within the disks of spiral galaxies. In particular, 
they find a dependence on stellar mass surface density, with 
the molecular fraction steadily increasing from surface mass 
densities of 10^'^ to 10^ Mq kpc~^. Our molecular fractions 
are smaller because our measurements are integrated over 
entire galaxies, but the same qualitative trend is observed 
for our global measurements. 

Even though fni and fH2 show different dependencies 
on /i* , the fraction of galaxies with non-detections as a func- 
tion of /i* in the HI and CO samples exhibi ts very similar 
behav iour. As shown in Figure [Sj) and in ICatinella et al.l 
(|2Q1QI 1. there is a critical mass surface density of /x* = 10 ' 
Mq kpc~^ below which all galaxies have a sizeable HI and 
H2 component, and above which cold gas seems to have 
mostly disappeared. There are similar thresholds in con- 
centration index and NUV— r colour (see e.g. our detec- 
tion fractions in Fig. [2). However, while the mean value of 
fHi ne ver falls below ^ 2 %, even at high /i* or NUV— r 
colour (!Fabello et alll2010h . /if2 drops sharply below that 
level, as shown in Fig. [71 and evidenced by the results of 
our stacking ex p erime nt (Fig. [6]). Furthermore, even though 
ICatinella et al.' (2010') found some red sequence galaxies 
with a surprisingly large HI component, none of these galax- 
ies have a sizeable molecular gas mass; none of the galaxies 
with NUV— r> 5 are securely detected in CO. It therefore 
seems that above our empirical thresholds in /i*, concen- 
tration index and NUV— r colour, an increasing fraction of 
galaxies with any form of cold gas at all appears to be dom- 
inated by atomic gas. This result is also striking in Fig. [9t, 
where we see the same population of atomic-gas dominated 
galaxies (i.e. the upper limits in Rmoi) almost exclusively 
at C > 2.6. We note that this critical value of concentra- 
tion index (C=2.6) corresponds to the ob served transition 



between late- and early-typ e galaxies (e.^. Shimasaku et al.l 
I2OOII : iNakamura et al.fi2003l : IWeinm ' 



et al.ll2009l l. 



Finally, Figure [9li shows that Rmoi is an increasing func- 
tion (r = 0.25) of N UV— r colour for galaxies with detections 
of both HI and CO. ICatinella et aL: (2010) reported a strong 
ant i- correlation between HI mass fraction and colour. Figure 
[SJl shows that there is also a fairly strong anti-correlation 
bewteen fH2 and NUV— r colour. The fact that the anti- 
correlation of /hi with NUV— r appears to "win" proba- 
bly reflects the fact that HI is dominant in regions such as 
outer galaxy disks where dust content is low and most of 
the starlight emitted by forming populations is emitted at 
UV wavelengths, whereas molecular gas tends to occur in 
the inner regions of galaxies where dust content is high and 
much of the light from young stars may be emitted in the 
infrared. In future work, we plan to look more carefully at 
these issues. 



7 COMPARISON WITH PREVIOUS WORK 

To put the new COLD CASS results in context, we assem- 
bled CO data from the literature for 263 nearby galaxies 
in the SPSS survey. Th e y are taken from t he compilations 
of iBettoni et al.1 (I2OO: ' 



- ,-■ ."T - 

('2004 h. lAlbrecht et al.1 (l2007l\. |Komugi et al." ("2008") and 



Yao et al 



(120031). [Casasola et al.1 



[K km s-^]-^) and 
In addition, the SDSS 



Obre schkow Rawlingsl (|2009l l . When multiple measure- 
ments are found for the same galaxy, the newest is assumed 
to supersede previous values. The values of were then 
homogenized to the best of our ability using a common con- 
version factor (Xco — 2.3 x 10^° cn 
cosmology {Hq — 70 km s~^ Mpc~^) 
photometry was reprocessed using the same technique used 
for the COLD CASS galaxies (see and used to measure 
reliable and homogeneous stellar masses. 

In Fig. llOb. we show the relation between and 
M* in this reference sample. The molecular gas mass frac- 
tion has a spread of more than two orders of magnitude, and 
does not appear to correlate with stellar mass. A significant 
fraction of this observed scatter can be attributed to mea- 
surement errors and inhomogeneities in the sample. The vast 
majority of these galaxies have z < 0.02, and therefore tend 
to be significantly larger than the typical observing beam, 
resulting in important aperture problems. Additional con- 
tributions to the artificially large scatter include different 
telescope calibrations, low S/N detections and selection on 
IR luminosity, which tends to bias fH2 high. 

As a comparison, we als o sho w as contours in Figure 
[TOk the relations produced bv lFu et al.l (|2010|) through semi- 
analytic modeling of galaxy formation including detailed 
prescriptions for the break-up of gas between the atomic 
and molecular phases. The models predict a significantly 
smaller range in Mh2/M^ than seen in the literature compi- 
lation. But it is also clear that a systematically measured set 
of galaxies will produce more consisten t results, for exam- 
ple the THINGS / HERACLES sample (|Walter et al.ll2008l : 
"Lero v et aDl2009h . which is also shown for comparison in 
Figure flQl 

The equivalent relation from COLD CASS is plotted 
in Figure llOb . showcasing the significantly reduced obser- 
vational scatter compared to the literature compilation, but 
the increased dynamic range due to the rigorously measured 
upper limits for the CO non-detections. We note that the 
COLD CASS and HERACLES galaxies span a similar re- 
gion in the plots. The HERACLES galaxies provide resolved 
CO maps for a much smaller sample of galaxies, so the two 
approaches are highly complementary. 

The main results presented in this paper are overall 
qua litativ ely similar to some earlier observations. For exam- 
ple. Sage (1993) found that MH2/Mdyn is independent of 
morphology, just like we find to be independent of 

concentration index. However, they found tha t Mh2/Mhi 
is a strong function of Hubble type (see also lYoung et al.l 
1989), which we do not. Based on these studies, we would 
have expected to find a significant population of early-type 
galaxies with large values of Mh2/Mhi- We see no evidence 
for such systems in our sample. 

Much of the earlier work regarding CO in nearby galax- 
ies focussed on the trends b etween Mh2, Mhi and m orphol- 



ogy, interaction state (e.g. iBraine Combeslll993l). or far 
infrared properties of the svstems (e.g. ISanders &: Mirabell 
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Figure 10. left: Relation between H2 mass fraction and stellar mass for a compilation of 263 galaxies with CO measurements 
available from the literature (gray filled circles). The values of have been homogenized for cosmology and the CO-to-if2 

conversion factor, but even with these correcti ons a scat t er of more than two orders of magnitude is observed. This is in strong 
contrast with pre dictions from the models of IFu etliD (|2010l ) (contours), and smaller but more homogeneous datasets (e.g. 
iLerov et al ] l2009l . shown as red circles), right: In comparison, we show the same plot for the COLD GASS sample, with our 
minimum stellar mass selection (vertical dashed line) and our integration limit (dotted line) indicated. 



ll985l V Using a set of galaxies still mostly based o n the 
infrared-based FCRAO sample, iBothwell etaP (|2009l l how- 
ever derive relations similar to ours, between gas fraction 
(atomic and molecular) and 5— band luminosity. Their find- 
ings are qualitatively similar to ours; they see that Mhi /M^ 
decreases with luminosity, but MuilM^ does not. The 
breakthrough here is that COLD GASS allows us to quanti- 
tatively describe, in an unbiased sample, how the molecular 
gas component varies with several key physical parameters 
which are at the basis of the theoretical effort towards un- 
derstanding the star formation process. 



8 SUMMARY 

We are conducting COLD GASS, a legacy survey for molec- 
ular gas in nearby galaxies. We target at least 350 massive 
galaxies (M* > 10^° M©) in the CO (1-0) emission line with 
the IRAM 30m telescope. Because the survey is unbiased, 
it will provide us with a complete view of the molecular gas 
properties of massive galaxies in the local universe, as well as 
the relations between molecular gas and other global galaxy 
properties. The stellar mass and redshift ranges also ensure 
that we recover the total CO line flux of the galaxies with a 
single pointing of the IRAM 30m telescope, and that a single 
CO luminosity to Mh2 conversion factor is likely adequate. 
Finally, our observations provide stringent upper limits on 
molecular gas fraction Mh2/M^ < 0.015 in the case of CO 
non-detections. 

In this paper, we present a catalog of CO (1-0) fluxes 
and H2 masses (or upper-limits) for the flrst 222 galaxies 
observed as part of COLD GASS, and report on their prop- 
erties: 



(i) The detected molecular gas mass fractions (Mh2/M^) 
are in the range of 0.9% to 20%, with a mean value < 
Mi/2/M* >= 0.066 ± 0.039. The mean gas mass fraction 
among the detected galaxies does not vary strongly with any 
global galaxy property except colour. 

(ii) We detect the CO (1-0) line in - 50% of COLD GASS 
galaxies, and while the detection rate is independent of stel- 
lar mass, it is a strongly decreasing function of stellar mass 
surface density, concentration index and NUV— r colour. 
None of the 68 galaxies redder than NUV— r=5 were de- 
tected, and stacking them leads to a non-detection and a 
stringent upper limit of < Mh2/M^ >= 0.0016 =b 0.0005. 

(iii) The mean molecular gas mass fraction (averaged over 
galaxies with detections and non-detections of the CO line) 
is a roughly constant function of stellar mass, but a de- 
creasing function of stellar mass surface density and con- 
centration i ndex. The observ ed trends are weaker than those 
observed bv lCatinella et all (2010) for the atomic gas mass 
fraction in a similar sample of galaxies. Of all parameters in- 
vestigated here, the molecular gas correlates most strongly 
with NUV— r colour, which is a tracer of speciflc star forma- 
tion. 

(iv) The molecular-to-atomic mass ratio, Rmoh hcis a 
mean value of 30% over the entire sample. It is a weakly 
increasing function of M* (r = 0.23) and /x* (r = 0.37). 

One result that we wish to highlight is the existence 
of sharp thresholds in galaxy structural parameters such 
as stellar surface mass density, below which most galaxies 
have measurable atomic and molecular gas components, but 
above which the detection rate of both the HI and CO lines 
drops dra stically. This result wa s discussed previously for 
the HI in ICatinella et al ] (|2010l ). The fact that the same 
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sharp thresholds also apply to the CO strongly suggests that 
onset of "quenching" processes in galaxies was associated 
with a change in their structure. We note that the same 
sharp drop in cold gas content is not seen as a function of 
stellar mass. Intriguingly, atomic gas dominates in the mi- 
nority of galaxies that are above threshold and that have 
significant cold gas content. One possible interpretation is 
that re-accretion of gas may still be possible following the 
quenching event. In future work, we will examine galaxies 
on either side of the "quenching threshold" in more detail. 
We will also look more closely at the relationhip between 
molecular gas content and star formation. 
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APPENDIX A: SPECTRAL GALLERY 

Al Galaxies with CO line detections 

In Figures lAlllASl we present SDSS imaging and IRAM 
spectra for the COLD GASS galaxies present in this data re- 
lease for which we have securely detected the CO (1-0) line. 
These are defined as those with S/N > 5 in the CO line. 
The tentative detections, those with S/N < 5, are shown 
in Figures IA9I and lAlOl All galaxies with a secure CO line 
detection have colours bluer that NUV— r=5.0, but cover 
uniformly the stellar mass range probed by COLD GASS. 

A2 Galaxies with non-detection of the CO line 

For completeness, we show in Figures IAllj|A13l the SDSS 
three-colour images of the CO non-detections. These objects 
tend to be red (NUV— r> 5), early-type looking (C > 2.6) 
galaxies. 

A3 Offset pointings 

As described in ^4.4[ we perform additional off-center point- 
ings for a small fraction of the COLD GASS sources. They 
are the larger galaxies {D25 > 40^^), with a strong detec- 
tion in the central pointing. These offsets are taken 16^^ (or 
three-quarters of the main beam) away from the centers. In 
Figures [Alil and [AT5l we show the SDSS images and IRAM 
spectra for the offset pointings. The data reduction process 
is identical to that used for the central pointings and de- 
scribed in Section 14.31 The emission lines are identified by 
examining the spectra at the expected positions based on 
the redshift on each galaxy, and the width of CO the line 
in the central pointing is used as additional information. 
Fluxes are measured by integrating over the region identi- 
fied though this process, and the flux ratio between offset 
and central pointings is used to determine an appropriate 
aperture correction for each galaxy using the technique de- 
scribed in Section [4.41 
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Figure Al. SDSS three colour image (1.5^ x 1.5^) and CO(l-O) line spectrum of COLD GASS targets with a secure detection 
(S/N > 5). The solid red line shows the systemic redshift of the galaxy as determined from the SDSS fiber spectra. The dashed 
blue line is the central velocity of the CO line, and the interval delimited by the two dotted blue lines is VKSOco, the full line 
width of the CO emission measured at half intensity. 
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Figure A2. continued from Figure IaT] 
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Figure A3, continued from Figure IaT] 
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Figure A4. continued from Figure IaT] 
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Figure A5. continued from Figure IaT] 
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Figure A6. continued from Figure IaT] 
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Figure A7. continued from Figure IaT] 
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Figure A9. SDSS three colour image and CO(l-O) line spectrum of COLD GASS targets with a tentative detection {S/N < 5). 
Lines are as described in the caption of Figure lAll Given the low S/N of these detections, the width determination is sometimes 
poor, and in these cases the corresponding vertical dotted lines are omitted. 
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Figure AlO. end from Figure |A9] 
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Figure All. SDSS three colour images of COLD GASS galaxies where the CO (1-0) line is not detected. Each image has a 
size of 1.5^ X 1.5^ 
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Figure A12. continued from Figure lAllI 
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Figure A13. end from Figure [Alll 
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Figure A14. SDSS three colour image and CO(l-O) line spectrum measured in offset pointings on COLD GASS galaxies. The 
blue circle shows the central position of the 22" beam, and the orange circle the offset position, located three quarters of a beam 
from the center (for G1977, G4216, G9551, Gil 112 and G42013, the offset was taken one full beam away from the center, see 
^4.4p . On the spectra, we overplot the optical redshift of the galaxies (red lines), and the central position and widths of the line 
detected in the central pointings (blue lines). The shaded region shows the region of the offset spectrum integrated to measure 
the line flux. The optical diameter from SDSS ^f— band imaging (1^25) and the flux ratio between offset and central pointings 
(/^jj) are given in each case. . 
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Figure A15. continued from Figure [A14I 
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